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The free energy, non-gradient terms of the Ginzburg-Landau (GL) expansion, and the 
jump of the specific heat of a multiband anisotropic-gap clean BCS superconductor are 
derived in the framework of a separable-kernel approximation. Results for a two-band 
superconductor, d-wave superconductor, and some recent models for MgB2 are worked 
out as special cases of the general approach. The classical results for the GL coefficients 
are derived in a simple way, directly from the general expression for the free energy of a 
BCS superconductor. 
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1. Introduction 



The Landau theory of second-order phase transitionai and its reahzation for 
superconductors, the Ginzburg-Landau (GL) gauge theory,a can be classified as be- 
longing to the most illuminating theoretical achievements in XXth-century physics. 
The basic concepts advanced in these theories often find applications in interdisci- 
plinary research. The microscopic Bardeen-Cooper-Schrieffer (BCS) theoryll makes 
it possible to calculate the parameters of the GL theory. Thus the phenomenology 
of superconductivity can be reliably derived once the parameters of the microscopic 
Hamiltonian are specified. Such a scheme ensures that there is no missing link 
between the microscopic theory and the material properties of the superconductors. 

The recent progress in physics of anisotropic-gap superconductors, e.g., cuprates, 
borocarbides and MgB2, has attracted considerable attention. However, the analysis 
of the thermodynamic behavior of these compounds appeared to be erroneous in 
some cases, and often the results from the classical papers in the field are not 
taken into account. Thus, the present situation created the necessity that some old 
theoretical results be delivered in a form suitable for fitting the experimental data. 

The purpose of the present paper is to provide a simple methodological deriva- 
tion of the non-gradient terms of the GL expansion. These are then used to set up 
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appropriate for numerical implementation expressions which employ the electron 
band dispersion and momentum dependence of the gap function at the Fermi sur- 
face. We have paid particular attention to reflect the relation of our results to those 
obtainedjn the pioneering papers in this field. Employing the Bogoliubov-Valatin 
methodcrO and the standard expression for the free energy in the HCS model we 
present a simple derivation of the results by Pokrovskii and RyvkinH and Gor'kov 
and Melik-BarkhudarovE We have particularly focused on the jump of the specific 
heat AC at the critical temperature Tc which is expressed through the coefficients 
of the GL expansion. Finally, the general formulae are applied to derive AC in 
different models used to describe the superconducting cuprates, borocarbides, and 
MgB2. 

2. Microscopic formalism 

Our starting point is the total reduced Hamiltoniani^i0 for a multiband super- 
conductor with singlet pairing 

H' = Ho-^iN + iJint = ^0 + ^int 

= ^''^P 4,paC6,pa + Jf'YYl ^b^P-y^P' 4,pT 4-piC6' .-p'iCfe' ,p'T ■ (1) 

b,p.a b,p b' .p' 

Here ^b,p = £b,p — /i is the quasi-particle spectrum of the "primed"^ non-interacting 
Hamiltonian Hq, with b being the band index, and p^^, Cb^pa are the quasi-particle 
creation and annihilation operators, respectively, for the &th band, quasimomentum 
p, and spin projection a =tii ■ The second, four-fermion term in Eq. (Q) is de- 
termined by the quasiparticle pairing interaction. The form of the corresponding 
pairing matrix elements Vf,^p;f,'^p' has been one of the vexing problems in the theory 
of high-Tc superconductivity as the underlying mechanism giving rise to pairing is 
still unknown. 

Howevfir^once the interaction is specified, one can apply the standard BCS 
treatmentall2l to obtain an equation for the superconducting gap. For the Hamilto- 
nian (|l]) this leads to the BCS gap equation of the following familiar form 

where 

Eb.P = ^^lp + \^b.p\^ nt^p = [exp{Eb^p/k,T) + 1]-^ (3) 

are the quasiparticle energies, and the Fermi filling factors, respectively, k-^ being 
the Boltzmann constant, and T the temperature. The summation over the band 
index b' is restricted to the partially filled (metallic) bands, containing the sheets 
which compose the Fermi surface. 

tAs usual it is convenient to work at fixed chemical potential fi and introduce Hq = Hq — fiN, 
where TV is the quasiparticle number operator. 
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Generally, the BCS gap equation requires numerical treatment, but for the spe- 
cial class of separable kernels V it can be reduced to a transcendental equation. This 
is the main reason why the gap symmetry of some exotic high-Tc superconductors is 
addressed on the basis of model separable kernels. For specific pairing mechanisms, 
however, it can be shownM that the interaction kernel is of separable form. As a 
rule simple model estimates result in a sum of such separable potentials. Therefore, 
as a next step we make the separability ansat£3 for the pairing interaction kernel 

Vfc,p;h',p' = -GXb.pXb'.p', (4) 

where G is a parameter, characteristic for the pairing interaction process, and the 
momentum dependence is determined by the (generally complex-valued) anisotropy 
function Xb,p- Due to the separable form of the BCS kernel (^) the temperature and 
momentum dependences of the gap also factorize, 

Ab,p(r) = S(T)x6,p. (5) 

As emphasized by Gor'kov and Melik-Barkhudarovl the factorization of the super- 
conducting gap into temperature (S) and momentum-dependent (x) multipliers is a 
general result due to Pokrovskii and Ryvkin,El valid for an arbitrary weak-coupling 
pairing interaction. The gap anisotropy function is given by the eigenfunction of 
the pairing kernel corresponding to the lowest eigenvalue. The factorable kernel 
Eq. (^) is just an interpolating Hamiltonian resulting in the same thermodynamic 
behavior as that determined by the initial kernel. Hence in what follows, without 
loss of generality, we shall use a separable interaction. Such a type of interaction, 
however, arises in lattice multiband models if one considers only local (single-site) 
interactions. 

Substituting Eq. (^ into Eq. (^) and passing from summation (for Af — s- oo) to 
integration according to the general rule in the Z?-dimensional case 

^E/(p) = f^-/"(^/(p)M/>p, (6) 

we obtain a simple equation for the temperature dependence of the gap, 

where the quasiparticlc spectrum in expanded form reads 

Et,p = [(e.,p - E^f + |S(r)p Ix^pH , (8) 

with Ep = ^ being the Fermi energy. For the ease of the following discussion we 
shall consider the simplest case of a single band and suppress the band index. This 
does not entail any restriction on the generality of the derivation. 
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3. Thermodynamic properties 

In order to implement the GL idea for representing the free energy as a function 
of the superconducting ordernarameter, -F(S), we shall employ the Bogoliubov- 
Valatin variational approach.atl Let us recall the main framework of this standard 
procedure. As a first step one carries out a transformation of the variables, intro- 
ducing new Fermi operators (for simplicity we will consider real kernels but the 
appropriate generalization can be easily performed) : 



pT (9) 
V-pi J V -sint^p cos6^p J \ c_px J 

Defining the BCS ground state 

|BCS) =n(cos0p + sm0pc^^clp^) |0), (10) 
p 

it is easily verified that ^p„ |BCS) ^ ^ Cp^ |0), and (BCS|BCS) = 1 = (0|0). 
The superconducting gap function is then used for a suitable parameterization of 
the Op angle, tan26'p = Ap/^p. Denoting Up = cos0p, Vp = sin0p, 

it is straightforward to verify that the condensation amplitude 2upi;p and coherence 
factor Up — Vp read, respectively, 

2upVp = sin 2(?p = ul~vl = cos 20p^^. (12) 

^p J^p 

Within the self-consistent BCS approximation the V'- "particles" can be considered 
as noninteracting, and their thermal-averaged number and entropy are given by the 
familiar expressions, respectively, 

np = {fipa) = (■(/'L'i/'pQ) = [exp{Ep/k^T) + 1]"\ (13) 



5(T) = -2fc3 V[nplnnp + (l-np)ln(l-np)]. (14) 



p 

Similarly, substituting in Eq. (|l|) 

) = ( «p -«p ^ ( ] (15) 

c-pi J V «p "p y V "^-p^ J 

we obtain for the expectation value of the reduced Hamiltonian with respect to the 
BCS ground state (|I^ 



{H')-2Y^^p[vl + {ul-vl)np]^^ 
p 



(16) 
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The minimization of the free energy 



F = {H') - TS 



(17) 



with respect to Ap leads to the gap equation (^. Then the substitution of the 
so obtained gap in Eq. ( [l^ ) gives the desired form of the minimal free energy. It 
should be noted that in order for us to derive the free energy as a function of the 
order parameter, we have to use Eq. (||) as ansatz, thereby considering the 

order parameter S as an independent variable for a fixed momentum dependence of 
the gap function Ap = Sxp. Thus, the desired function _F(S,r) is obtained from 
Eq. by substituting the expression for the averaged energy, Eq. (|l^), and that 
for the entropy, Eq. (M), which can be also rewritten in the form 



Er^ tanh ■ ^ 



2k T 



2KTln 2 cosh 



2k T 



(18) 



The self-consistent BCS approximation gives an analytical dependence on the order 
parameter Ap — Exp, 



= 2 ^, 




1 



py Ep cxp{Ep/k^T) + l 

2 



tanh 


Ep 






2k,T 


>, 




Ep 


tanh 


Ep 


2k^T 


2k 


T 

B ^ 



In 2 cosh ■ 



Er. 



2k T 



(19) 



Then the Taylor expansion provides the coefficients in the GL functional per unit 
cell 



which is often a satisfactory approximation even far from the critical temperature. 
Let us mention that close to Tc 



a{T) w ao 



T-T, 



b{T) « b = b{Tc) > 0. 



(21) 



The order parameter S, being complex in the general case, is proportional to the 
GL effective wave function, while 77 is proportional to the superfluid density. 

A straightforward way to obtain the ao and b coefficients is to compare the 
variations with respect to rj of the expression ( |l7| ) and the GL expansion (|20|). 
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Introducing the notation f3 — l/k^T the former reads 

SF /eM' \ 



Srj 



-G 



IXpI 



■(1 - 2np) 



(^(l-2.p) 



-277 



IXpT 



[l-2np-2/3Spnp(l-np)] 



pj 



Furthermore, introducing the functions 



IXpl 

2Er 



.(l-2np; 



e(77, T) = ( ^ [(1 - 2np)e^ + 2/3£;p|App np(l - n^)] 



2EI 



Eq. ( p2D can be cast in the compact form 

5F 



Sr] 



B{ri,T){l-GA{ii,T)). 



(22) 

(23) 
(24) 

(25) 



Notice that the extremum condition of ([25|), 1— G^(r/, T) = 0, gives the gap equation 
(^ . Now varying Eq. ( pO| ) with respect to ?7 and comparing with Eq. (|2^) we find 



9-4(77, T) 



9r 



b = - 



dA{v,T) 



ri=0, T=T^ 



(26) 



,,=0, T=Tc 



where we have used the identity 6(0, Tc) = -4(0, Tc) = 1/G. Thus, Eq. (gO|) with 
and 



Eqs. (|2l|), (23), (|24D, and ( |26| ) provide the complete set of equations which deter- 
mines the GL free energy of an anisotropic-gap superconductor. 

We shall proceed now with working out explicit expressions on the basis of the 
general relations (p6|). Taking the corresponding derivatives of A{r],T) the oq and 
b coefficients read 



1 / IXpP 
Ak^Tc \cosh^ fp 



where Vp = S^p/2k^Tc and (cf. Ref. |K 



^, , if tanha; 



1 



1 



9 \ 

X'^ \ X 



cosh X 



,^,^r.f{\x.\'Q^^.))p 



tanh^ X tanhx — x 



(27) 



(28) 



In calculating the momentum averages (|22) one has to take into account that the 
integrands in the momentum-space integrals exhibit sharp maxima at the Fermi sur- 
face Ep = Ep. Thus, within acceptable accuracy, one can carry out an integration 
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along the normal to the Fermi surface whereupon fp may take on values in the inter- 
val (—00, +00) while the longitudinal quasimomentum component p|| varies negligi- 
bly. Taking into account the numerical values of the integrals {cosh x)^^ dx = 2, 
and 



16^ 1 14^^3^^ 



where C(™) is the Riemann zeta function,"'' we obtain the desired form 



(29) 



""^(i^'O.^ ''^8;;^SfvF(ixpa^ (30) 

Equation ( ^ ) reproduces the results by Pokrovskii and Ryvkini and Gor'kov and 
Melik-Barkhudarov0 obtained for the case of a clean superconductor. The detailed 
consideration of the influence of nonmagnetic impuiities on the anisotropicLXjrder 
parameter was given by Pokrovsky and PokrovskyJlil Very recently KoganEj has 
analyzed the macroscopic anisotropy in superconductors with anisotropic gap. 

For the momentum-space integrals over the Fermi surface which appear in 
Eq. (pO|) we use the notation 



(/(P)). . E / " ■ 7 - ^ri/(p)^ -111- 1 ^"'^ 



(31) 



where dSb,p is an infinitesimal surface element of the (possibly fragmented) Fermi 
surface sheet of the 6th energy band, and Vf,.p — VpEb^p is the quasiparticlc "veloc- 
ity" which according to the present convention has dimension of energy. Conversion 
to the true velocity can be performed by multiplying the (dimensionless) quasimo- 
mentum by h/a, where a is the lattice constant. For a conventional electron-phonon 
pairing mechanism the isotropic gap, Xp = I7 is often a reasonable approximation. 
Thereby averaging in Eq. ( ^0[ ) results in the density of states (per spin orientation 
a) at the Fermi level pF = (l)p ■ In this importaiit case the general formulae are 
in agreement with the classical result by Gor'kovlij (in particular, cf. Eq. (11) in 
Ref. 

In thermodynamics of second-order phase transitioiis the ratio of the GL coeffi- 
cients flQ and b determines the jump of the specific heatEl per unit cell at the critical 
temperature, so for superconductors 



tThis result is easily obtained by employing the infinite-series representation tanh(x/2) = 
4a; J]^Q [7r-^(2n -|- 1)^ -|- a;-^] ^ , and recalling the relation between the Hurwitz and the Riemann 
zeta functions, respectively, |) = EJT^o + 3)™ = (2™ " 1)C(™-)- 
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where Cg is the specific heat of the superconducting phase, and 



^l^^^klpFT, (33) 



is the normal-phase specific heat per unit cell at Tc; the factor 2 takes into account 
the spin degeneracy of the normal paramagnetic phase. These formulae can be 
derived directly from the BCS expression for the specific heatcfEU 

C(T) ^dS{T) 2 / f , T d ,\\ 

Substituting Eq. ( pO[ ) into ( |3^ ) , and taking into account that the order parameter 
S is not affected by momentum averaging, one finds the general expression for the 
relative jump of the specific heat of an anisotropic-gap superconductor 



AC 

Cn 



7C(3) /3a ' (| Apl^), ^ ^' 7C(3) " '-^'''^ ' ' ' ' ^^'^ 



the latter being the universal BCS ratio. § Note that /3a is similar to the Abrikosov 
parametenla /3a. 

4. Applications 



4.1. Jump of the specific heat for layered cuprates 

In models for layered cuprates one often postulates the following functional form 
of the gap anisotropy 

Ap cx cos(2(/)p), 

where (j)p = arctan(pj,/pa;). In the case of a parabolic dispersion, oc p^, averaging 
in Eq. (j3|) is straightforward and yields 

1//3a = 2/3. 

This value should be considered as a lower bound for 1//3a as the Fermi velocity 
is minimal where the supejmonducting gap is maximal. Consistent with that, for 
a realistic band dispersionll3 and gap anisotropyllil the reduction factor may well 
reach 1//3a = 0.72. 

As pointed out by AbrikosovJ^^ if the Fermi level is very close to an extended 
saddle-point singularity the derivation of the GL equations requires a special con- 
sideration. 

§The numerical value for the specific heat "anomaly" at Tc cited originally by Bardeen, Cooper 
and SchriefTerg is slightly different, AC/Cnh, = 1-52. 
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Figure 1: Jump of the specific heat AC/Cn\Tc versus the "equatorial-to-polar" gap ratio 
Ae/Ap. For a given AC/Cn|p^ value we have oblate Ae/Ap > 1 and prolate A^/Ap < 
1 solutions. The solid curvecj represents Eq. (|36|), and is applicable to the model 1^ 
Posazhennikova, Dahm and MakiEII The dash-dotted, line is our analytical solution,Ej 
Eq. (^tI), applicable for the model by Haas and Maki.tj The dashed liiteiis the jump ratio 
AC/Cn|Tc — 0.82 ± 10% measured by Wang, Plackowski, and Junod,E3 with the shaded 
area showing the experimental error bar. 




4.2. Jump of the specific heat for MgB2 



Let us now make some eomparison between different models applied to MgB2. 
In a very recent e-print Posazhennikova, Dahm and Makicil discuss a model for 
the gap anisotropy in MgB2, a material which has attracted a lot of attention 
from «3ndensed-matter physicists in the past two years. A central issue in this 
workEil is to propose an analytic model for analyzing thermodynamic behavior. 
Assuming a spherical Fermi surface, a simple gap anisotropy function is suggested, 
A(p) ~ Ag/Vl + Az^, where z = cos9, and 9 is the polar angle. This model 
leads to useful results for the temperature dependence of the upper critical field 
Hc2 and of the specific heat, which can be fitted to the experimental data, thereby 
determining the optimal anisotropy parameter A. Note that A — (Ag/Ap)^ — 1, 
with Ap = A(z = 1) and Ag = A(z = 0), and the gap ratio is parameterized as 
Ae/Ap = VTTA > 0. 

We shall now apply the general results obtained in Sec. 3 to derive a conve- 
nient analytical expression giving the possibility for determining Ag/Ap from the 
available data for the jump of the specific heat.E3 Following the weak-coupling BCS 
approach j£ilo from Eq. we fincO for ^ > 0, and — 1 < A < 0, respectively, 



ACiA) 



12 2(l + w^)(arctanw)^ ,„ - ./T - 
i(,[o) w + w(l + w ) arctanw 



12 2{l~y^){tanh-'yy 



2/(1 - y ) tanh y' 



y = iw — \/—A = \ 1 — 



A. 



(36) 
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For a given specific-heat jump, this expression leads to two solutions (oblate, Ag/Ap > 
1, and prolate, Ag/Ap < 1). The relevant example is shown in Fig. |^. In Ref. ^ 
values of AC{A)/Cn\Ta consislept with Eq. ( ^6| ) have been tabulated. The anal- 
ysis of the experimental datacScZl for the angular dependence of i/c2 unambigu- 
ously demonstrates that one has to analyze only the "oblate" case. The exper- 
imental value AC/Cnlr, = 0.82 ± 10%, reported in Ref. E3, gives ^ ~ 16 and 
Ae/Ap « w « 4.0 ± 10%. For this significant anisotropy, the "distribution" of 
Cooper pairs (|Ap| ) oc l/[p1 + {pf/w)'^] has a maximum at — 0, and pp is 
the Fermi momentum. This general qualitative conclusion is in agreement with the 
hints from band calculations that the maximal order parameter is concentrated in 
a nearly two-dimensional electron band, but all bands Sb^p contribute to Cn/M- 

In another paper Haas and MakiEj coxLsidered the model gap anisotropy A (x 
1 + az^ for which similar calculation giveaH3 



AC 



12 l + 4a/3 + 38aV45 + 4a3/15 + 0^/25 
7C(3) l-H4a/3-f 6a2/5-H4aV7-l-aV9 ' 



This model, however, cannot explain the significant reduction of AC/ChItc observed 
experimentally, as can be seen in Fig. 0. 

4.2.1. AC/CiiItc within a two-band model 

For the Moskalenkoi-two-band model, advocated for the first time for MgB2 
in Ref. Eq, Jilq. (p5|) gives (to within a typographical correction) the result by 



MoskalenkoJa 



AC 



where 



12 (|Ai|Vi + |A2pp2)^ 
7C(3) (pi+p2)(|Ai|V + |A2|V2) 

. 1.426l5i±ii^< 1.43, (38) 



A2 Pl + P2 

and Pl and p2 are the densities of states for the two bands at the Fermi level. 
In Fig. |2| we present the specific heat jump as a function of 6 for various values 
of X. Taking for an illustration x = 0.515 and AC/ChIto = 0.82, Eq. ( ^ gives 
A1/A2 « 4.0 in agreement with the maximal-to-minimal-gap ratio Ag/Ap « 4.0 
obtained using Eq. (^. 

We should stress here that the two-band model also describes a kind of "anisotro- 
py" in the sense that a ^^non-constanV^ order parameter, having different constant 
values within the different bands comprising the Fermi surface, leads to modified 
GL parameters and reduces the jump of the specific heat. For a survey on a set 
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Figure 2: Jump of the specific heat AC/Cn at Tc versus the gap ratio A1/A2 for a set of 
different x values. The latter parameter is the relative contribution of the "first" band to 
the density of states within the two-band model for MgB2. 



of parameters see Table I in Ref. g9[ Certainly, from the jump of the specific heat 
alone one cannot judge about the validity of any model, so subtleties, e.g., related 
to strong coupling effects and other anisotropics, can be hidden in the spread of the 
parameters in the table mentioned. 

As we have aimed here on methodological aspects as well, it is beyond the scope 
of this work to analyze in detail different experimental data and the theoretical fits 
to them. Nonetheless we shall mention that the Moskalenko-two-band modelcl with 
isotropic gaps agrees betteiEj with the experimental data for MgB2 compared to 
other proposed models, but it would be premature to make any final judgment. 
Thus the synthesis of crystals with really sharp transition may well be rendered 
necessary. 



4.3. Effect of finite Debye frequency tJD for phonon superconductors 

Up to now we have dealt with energy scales corresponding to very high phonon 
frequencies or exchange interaction. For phonon superconductors, however, we have 
to take into account only a narrow "layer" of pairing electrons with energies |^p| < 
fkuY). In this case Up G {—huj-a /2k^Tc, +fkj-£,/2kgTc), and according to Eq. (|30| ) we 
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have 



F 2k^Tc' 



b-^TlA^ilXpf)^ I Qi^)dx. (39) 



'8ik,T,y /f 

For the relative jump of the specific heat, using Eq. (p^) and Eq. (p5[), we obtain 



AC 



12 1 ta.nh^{hiOD/2k^Tc) 



(40) 



Kishore and Lamball3 have recently shown that the behavior of AC'/Tc calculated 
wilhin the BCS model with finite uu is very similar to the results by Marsiglio et 
a/B3 and CarbotteE^I obtained from the Eliashberg theory; see also Ref. p2_jand 
Refs. As regards the recent first-principles linear-response calculationla for 

MgB2 the agreement between the theory and experiment is comparable to the qual- 



ity of the fitcj to experimental data employing the iiriginal week-coupling two-band 
model.u Precise analysis, however, has manifestedES a significant temperature de- 
pendence of A1/A2, which we interpret as an evidence for strong coupling effects. 
Further support can be possibly gained from the direct comparison between the 
calculated relative jump and the week coupling formula Eq. ^ 

5. Summary 

Analyzing the Ginzburg-Landau coefficients and the jump of the specific heat 
at Tc we were able to establish that a number of special cases discussed in the 
recent literature on superconducting cuprates, MgB2 and borocarbides can be easily 
derived as a consequence of the classical treatment of thermodynamics of anisotropic 
superconductors carried out in the dawn of the BCS era. It is somewhat strange that 
all those early papers are not being used in the current literature and a number of 
results are rediscovered or treated by a numerical brute force. We have cast, where 
necessary with the proper correction(s), those classical results into a form to be 
used for fitting the experimental data. The relative jump of the specific heat is one 
of the most appropriate quantities which can provide important information for the 
properties of superconductors. Thus, it wouki be interesting to see an analogous 
approach for unconventional superconductorsE3 as well. 
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